Measuring Spin Accumulations with Current Noise 
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We investigate the time-dependent fluctuations of the electric current injected from a reservoir 
with a non-equilibrium spin accumulation into a mesoscopic conductor. We show how the current 
noise power directly reflects the magnitude of the spin accumulation in two easily noticeable ways. 
First, as the temperature is lowered, the small-bias noise saturates at a value determined by the spin 
accumulation. Second, in the presence of spin-orbit interactions in the conductor, the current noise 
exhibits a sample-dependent mesoscopic asymmetry under reversal of the electric current direction. 
These features provide for a purely electric protocol for measuring spin accumulations. 



PACS numbers: 73.23.-b, 72.25.Dc, 85.75.-d 

Noise measurements on non-equilibrium electric cur- 
rents are very efficient probes of the dynamics and na- 
ture of the charge carriers [T]. At low temperature, the 
classical Johnson-Nyquist noise is suppressed and quan- 
tum effects govern the behavior of the surviving shot 
noise. In the mesoscopic regime, the noise power 5" is re- 
duced below its uncorrelated Poisson value 5*0 = 2\q\{I), 
where (/) is the average electric current, by the Fano 
factor F = S/Sq. The value of F depends on the elec- 
tronic dynamics. For instance, one finds _F = 1/3 in 
diffusive systems and F = 1/4 in ballistic chaotic sys- 
tems pQ [5] . Alternatively, shot noise measurements have 
determined the charge \q\ of current-carrying quasiparti- 
cles in normal-metal/superconductor junctions and in the 
fractional quantum Hall effect ^{3]. In this manuscript 
we further illustrate the usefulness of current noise mea- 
surements by showing how they can reveal the magni- 
tude of non-equilibrium spin accumulations. Our results 
provide for a purely electric protocol to measure spin ac- 
cumulations, which has the potential to quantitatively 
determine their magnitude. It therefore goes one step 
further than the optical methods used so far to detect 
magneto-electrically generated spin accumulations [H [S] . 
Alternatively, the noise measurement we propose, cou- 
pled with an electric measurement of the spin Hall and 
inverse spin Hall effects [BHH] , can provide key experimen- 
tal information on the conversion between spin accumu- 
lations and spin currents. 

A number of works have investigated charge current 
noise from polarized reservoirs. Reference [3] suggested 
using current and noise measurements in the single- 
channel limit to measure the spin injection efficiency from 
a ferromagnet for weak spin flip scattering. Other re- 
lated works have pointed out that noise measurements 
in hybrid paramagnetic/ferromagnetic structures can re- 
veal information on the relative orientation of the fcrro- 
magnets [lOj and on the spin relaxation processes in the 
paramagnet pTMl4j . These results have been at least 
partially confirmed by numerical simulations |15) . In 




FIG. 1: Three-terminal quantum dot connected to two 
unpolarized electron reservoirs (labeled 1 and 2) and one 
reservoir (3) with a non-equilibrium spin accumulation. Top 
right inset: a ferromagnetic reservoir with an equilibrium 
spin accumulation (a case we do not consider in this paper). 

non-interacting systems, current cross-correlations have 
a sign determined by the statistics of the charge carri- 
ers. Investigations of a single-level interacting fermionic 
quantum dot coupled to ferromagnetic leads have demon- 
strated the emergence of positive (boson like) current 
cross-correlations for certain relative orientations of the 
polarizations [16) . In all these instances, only ferromag- 
netic, i.e., equilibrium polarizations were considered. 
Below we show that non-equilibrium spin accumulations 
generate fundamentally different electric current noises. 
Our main findings are that (i) at low enough tempera- 
ture, the small-bias noise saturates at a value reflecting 
the spin accumulation, and (ii) in the presence of spin- 
orbit interactions, the current noise exhibits a sample- 
dependent, mesoscopic asymmetry under reversal of the 
electric current direction. These two features appear only 
in the presence of non-equilibrium spin accumulations. 

We consider a system such as the one sketched in 
Fig. [1} where a mesoscopic conductor is connected 
via multichannel leads to M external reservoirs, a — 
1,2,..., M, at electro-chemical potentials /Xq, = {fia^ + 
/iQ^)/2 and with non-equilibrium spin accumulations 
^fJ-a — (MqT ^/^a^)/2, along reservoir-dependent axes de- 
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fined by unit vectors uia = imax-,may,iTT'az)- We use the finear response scattering approach to transport to write 
the zero-frequency noise power in units of (e^/h) as [II I15j 
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where is the Fermi function for electrons with spin cr = ± along in terminal 7, and the sums run over all 
terminals 7 and 5 (including a and all channels m G 7 and n £ 5, and all spin orientations a, a' = ±. We defined 



m,cr:n,cr 7 
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where Sa-y denotes the 2Na x 2N^ subblock of the scattering matrix of the total system, corresponding to scattering 
from lead 7 to lead a, N^^j being the number of channels in those leads. This assumes that is spin-independent 
in all leads, and we will comment on the case N^^ 7^ N^^ later. Equation ([T]) differs from Eq. (52) in Ref. Ij in that 
spin indices are explicitly written down here. All our calculations below are current-conserving, gauge invariant, and 
satisfy linear response reciprocity relations, as they should. 



We assume that the temperature, applied voltages, and 
spin accumulations are low enough that the scattering 
matrix is essentially constant in the energy interval where 
the square bracket in Eq. ([!]) does not vanish. We then 
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(1 - ff) + ff (1 - /-) 

and introduce the two-terminal symmetry coefficients 

-cSS ^ \^ -pcrcr' -pAA ^ \^ -paa' lA,,\ 

■^is = 4 -^-rS : J'-rS ^ll^'^'^ -^is , (4a) 

aa' aa' 

- \ E --^7^' , -^7"/ = \ E -'^'J ■ (4b) 

era' aa' 

The indices S (A) indicate that the function is symmet- 
ric (antisymmetric) with respect to the spin accumula- 
tion in the corresponding lead, e.g. J^^g^{5^^,5^s) — 
J^^f{-6fj,^,Sfis) = -J^^s^{Sfj,j,-6fis)- We obtain 



Safs = 2A;bT 



(5) 
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with the spin-dependent noise coefficients 



T-ai 



^ Tr 
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Here, the trace runs over both spin and channel indices, 
1^ is the Nj X identity matrix, = cr ■ m^, where cr 
is the vector of Pauli matrices, and cr^ is the 2x2 iden- 
tity matrix. The coefficients given by Eq. ^ generalize 
those introduced in Ref. [T7] for the calculation of spin 
conductance, to the calculation of noise. The linear re- 
sponse Eq. ([5]) is valid for any number of terminals whose 
temperatures, electro-chemical potentials, and spin accu- 
mulations are encoded in the coefficients and for any 
particle dynamics contained in the noise coefficients T. 



We first mention symmetry properties of the coeffi- 
cients T. Aside from their symmetry with respect to 
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spin accumulations [see Eqs. (|4|], they satisfy (i) T^"^ 



T-AA _ -pAS 
-^75 — -^(57 



if = 0, (u) J-f/ = J-//, (iu) J-ff 
and T^g^ are symmetric, while F^f and are an- 

tisymmetric with respect to the voltage bias between 7 
and (5, and (iv) = 0. Property (in) is of particular 
interest, since together with Eq. ([s]), it implies that in the 
presence of spin-orbit interactions, the noise power is no 
longer symmetric under reversal of the current/voltage 
when there is spin accumulation in at least one reservoir. 

The system-dependent noise coefficients T are deter- 
mined by the orbital and spin dynamics of the electrons. 
We calculate their mesoscopic ensemble average and, 
when it vanishes, their typical value, taken as the root 
mean square of their distribution. In the absence of spin 
accumulation, only spin-independent coefficients T™^^ 
enter Eq. ([5]), whose mesoscopic averages (T™^^) have 
been computed using, for example, random matrix the- 
ory [TH| or the trajectory-based semiclassical theory . 
Extended to account for the Pauli matrices in Eq. (|6|, 
these methods give for chaotic ballistic systems 
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a result which holds to leading order in the total num- 
ber of channels, A^t = J2a^a ^ 1' both the 
unitary (broken time reversal symmetry) and the sym- 
plectic (broken spin rotational symmetry but preserved 
time reversal symmetry) ensembles |20| . In the orthogo- 
nal ensemble (preserved spin rotational and time rever- 
sal symmetries), Eq. Q holds provided one substitutes 
SaO ~^ 1- Then, in the case of non-collinear spin accumu- 
lations in leads 7 and 6, the symbol Sab for a = z — b 
should be understood as m.^ • mg. 

As a first example, we consider a spin preserving sys- 
tem with only collinear spin accumulations. This gives 
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75a/3 - T°Sai3 ^nd T^SaB = 0' O^ly spin-diagonal co- 
efficients enter Eq. (l5| and the two spin species are 
uncorrelated, with, additive contributions to the current 
noise. Despite zero charge current, the current noise can 
be finite in the presence of spin accumulations. 

Aiming at an all-electrical measurement protocol for 
spin accumulations, we show how the previous a pri- 
ori trivial observation carries over to spin systems with 
fully broken spin rotational symmetry, where the elec- 
tron dwell time is larger than the spin-orbit time. For 
simplicity, we focus on symmetric two-terminal geome- 
tries, N — Nl — Nji, with a spin accumulation only in 
the left lead, Sfii^ = Sfi 7^ 0, Sfi^ = 0, and with an applied 
voltage eV = /iL^MR- Current conservation ensures that 
S-RR = ^LL = -<5'rl = -S'lr, and accordingly we only 
discuss S = S'rr from now on. Equation ^ gives 

•^LR = (eV^ + crSfi) coth [(eV^ + aS^i) /2kBT] , (8a) 

= {a-a') 5 fi coth [{a ~ a') 5 fi/2kBT] , (8b) 

-^RL - J'Ik: J-^S' =2fcBr, (8c) 
while Eq. ([t]) gives 

('^7(5Rr) chaotic 

= {N/4)dab {Sao + 26^s) ■ (9) 

In the limit of zero temperature, we get the ensemble 
averaged zero-frequency noise power as 

(S) =(1/4)7V {\eV + 5fi\ + \eV - 5^1 + I-^mI) ■ (10) 

This function is plotted in Fig. 2(a). The spin accumu- 
lation manifests itself as a change in the slope of the 
noise at a crossover voltage \eV\ = \Sfi\, with a satura- 
tion at (S) = 3/4 X N\Sn\ for \eV\ < \Sn\, turning into 
(5) = A^/4x(2|ey| + |(5/i|) for \eV\ > \Sfi\. For Sfi = 0, we 
reproduce the result (S) — 2eFI with F = 1/4, valid for 
chaotic ballistic systems IJ. The abrupt change in slope 
at \eV\ = |(5/i| is smoothed out at finite temperature. 
This is shown in Fig. 2(b), where we plot the finite tem- 
perature analytic formula for (S) obtained from Eqs. ([5]), 
([S]), and (|9|. The crossover from low bias, \eV\ < \S^\, to 
high bias, \eV\ > \S^\, is still extractable from d^S/dV^, 
as illustrated in Fig. 2(c). The second derivative reaches 
its maximum close to \eV\ — \Sfi\ as long as k-^T < 
For zero applied voltage, F = 0, we get 



{S) ^{11/ A)NkBT + {l/A)N5ii 

X [2coth((5^/2/cBT) +coth((5^/fcBT)] 



(11) 



In the low temperature limit, k^T <^ \SfJ-\, the noise 
due to the spin accumulation decouples from the thermal 
noise, allowing for the measurement oi Sfj,hy varying the 
temperature. In the opposite limit, k^T 3> \Sfj,\, we re- 
cover the standard result for the Johnson-Nyquist noise, 
S = AkBTG. 

So far we have shown how a spin accumulation can 
be quantitatively extracted from the ensemble averaged 
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FIG. 2; Current noise in a two terminal conductor vs. applied 
bias voltage for a spin accumulation of 5/1 = 400 /leV in a 
single lead (solid lines) or no spin accumulation in either lead 
(dashed lines), (a) T = K and iV = 10. (b) T = 1 K 
and A'' — 10. (c) Second derivative of the data in panel (b). 
(d) Typical asymmetry in the current noise as a function of 
applied bias for T — 1 K and N — 2. 



current noise. According to Eq. ([t]), the average (T^ppirj) 
vanishes. However, individual samples might exhibit a 
nonzero T£|^j^j^, which, quite importantly, generates a 
contribution to the noise that is antisymmetric in the 
bias voltage. Using Eq. ^ we get, at zero temperature, 



SS ^ SiV) - Si-V) = 2 r£ V{ ley + 5fi\ - \eV 



(12) 



while at high temperature the effect is washed out, as 
expected: 6SiV) = 4/3 x T^^^r x 6neV/kBT. 

We estimate the magnitude of this asymmetry in a 
typical mesoscopic sample by calculating the root mean 
square of T^^^p- Again, using the method of Ref. [T5] . 
we find that in chaotic ballistic systems 



(varr£i^RR)chaotic = 1/128 + 0{N^^) 



(13) 



Accordingly, one has a typical asymmetry of 5Styp ~ 
Tms{6S) — \eV\/2^/2 at low voltages and SStyp ~ 
|5/i|/2-\/2 at higher voltages. This typical noise asym- 
metry is illustrated in Fig. 2(d). Interestingly, the asym- 
metry renders the noise smaller at finite voltage than at 
V = 0. A noise asymmetry was reported in Ref. [H] 
in systems with broken time-reversal symmetry in the 
nonlinear regime. The mechanism for this asymmetry is, 
however, different here. 

Because the asymmetry does not scale with the number 
of channels, while the total noise does, we predict that it 
is more evident in systems with few channels. The next 
order contributions tend to somewhat reduce the lead- 
ing order result in Eq. (13 1. This is most pronounced 



at iV = 1, where time-reversal symmetry requires that 
"^LRRR vanish identically. This is analogous to the van- 



ishing of Tr 



/3 « ) "^ipsap 



found in Ref. [2T]. Our 
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calculations therefore suggest that the asymmetry is best 
visible for TV = 2. 

The method of Ref. [18] can also be applied to diffusive 
systems with an elastic mean free path much smaller than 
the linear system size, £ L. One obtains 



dephasing, the noise asymmetry determined by Eq. ( 13 ) 



('^LLRr) diffusive 

/TOO \ 

\ ' LRRR/ diffusive 

(var Tlrrr ) diffusive 



{2/3)Ne/L, 
(2/35)£/L. 



(14) 
(15) 
(16) 



This gives, in particular, for T — 



(S) = {2Nl/iL) {\eV + VI + \eV - 6^i\ + 2|(5/i|) , (17) 



and for y = 

(5) ={ANl/ZL){?,kYiT 



X [coth((5^/2fcBT) + coth((5^/fcBT)]} . 



(18) 



Comparing Eqs. (10) and (11) with Eqs. (17) and (18) 



we see that after the substitution N — >■ N£/L, the noise 
averages for chaotic and diffusive conductors differ only 
by prefactors of order one. 

With the above results, we now evaluate the ratio of a 
typical noise asymmetry to the ensemble averaged noise. 
At \eV\ = \5fi\, where this ratio is maximal, we get, at 
zero temperature. 



SStyp/{S) = (l/N) X 



72/3, 
y/9L/70i, 



chaotic , 
diffusive . 



(19) 



Because metallic diffusive wires have TV ^ L/£, we see 
that a chaotic system is better suited for detection of spin 
accumulation from the noise asymmetry. 

We finally comment on the case of a spin dependent 
number of channels, Na^ ^ N^i, which occurs for large 
enough spin accumulations, 5pLa/^ia > ^/Na and breaks 
time-reversal symmetry. Equation ([7| becomes 



j^aj^b ^ _ 
N^il-Sab)]} 



(20) 



with TV. 



O/z 
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Interestingly, Eq. ( pO| ) implies 
a finite average asymmetry (SS) = 0{N^). 

In our derivation of Eq. ([s]), we neglected the energy 
dependence of the scattering matrix. This is legitimate 
as long as the expression in brackets in Eq. ([I]) is finite 
only in a narrow energy range. When this is not the case, 
the noise asymmetry will be damped even in individual 
samples SS typ — 0, unless the spin accumulation is large 
enough that TVl-]- ^ TVl^. Simultaneously, Eq. ^ may 
still give (S) provided one substitutes T — ^ (T). This is 
legitimate as long as the response is linear, meaning the 
applied voltages do not change the electrostatic profile of 
the conductor, and no substantial energy relaxation takes 
place in the system. We finally note that, in presence of 



is algebraically damped, in the same way as conductance 
fiuctuations are. [23] We thus believe that the noise asym- 
metry we predict is observable even when dephasing is 
taken into account. 
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